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Comment on “Accelerating cosmological expan¬ 
sion from shear and bulk viscosity” 

In a recent Letter [l| the cause of the acceleration of 
the present Universe has been identified with the shear 
viscosity of an imperfect relativistic fluid even in the ab¬ 
sence of any bulk viscous contribution. If C and rj denote, 
respectively, the bulk and the shear viscosity coefficients, 
the conclusion of Ref. [ij is that the deceleration param¬ 
eter can be negative when C 0 and rj ^ 0. The gist 
of this comment is that the shear viscosity, if anything, 
can only lead to an accelerated expansion over sufficiently 
small scales well inside the Hubble radius. 

When the fluctuations of the metric are perturbative, 
as in Ref. [ij , the gauge of the spatial gradients is given, 
in practice, by the product of the wavenumber k and 
of the conformal time coordinate r ~ l/'H. Since kr 
roughly equals the ratio between the Hubble radius and 
the typical spatial gradient. Ref. [Ij assumes that the 
Newtonian potentials are 0{'H?/k'^), the velocity field is 
0{'H/k) while the density (and its contrast) are both 
0{1). We remark that this approximation scheme is ten¬ 
able provided the relevant Fourier modes are located suf¬ 
ficiently inside the Hubble radius (i.e. kr I), as sug¬ 
gested in the previous paragraph. From this hierarchy all 
the equations easily follow if we set to zero the shift vec¬ 
tor in the four-dimensional space-time metric 0. On a 
more technical ground defining , the condi¬ 

tion gij^^u^u" = 1 for the four-velocity of the fluid implies 
u° = vT -\- /N which calls for an expansion in powers 
of u^, for instance, in the explicit expressions of the co¬ 
variant conservation equations [ij. 

The simplest way to proceed is to drop the assump¬ 
tion that the gravitational and hydrodynamic fluctua¬ 
tions are perturbatively small and study the effects of 
the bulk and shear viscosity within the expansion in 
spatial gradients. This nonperturbative scheme for the 
dynamics of the geometry and of the sources is the 
gravitational analog of companion methods employed in 
quantum field theory, in electrodynamics and in plasma 
physics 0. Order by order in the spatial gradients, the 
validity of the momentum and of the Hamiltonian con¬ 
straints can be explicitly enforced. In the synchronous 
frame 0 the nonlinear generalization of the decelera¬ 
tion parameter is given by q{t, x) = —1 + K/TrK^ where 
Kij = —^ijl2 is the extrinsic curvature and the overdot 
denotes a derivation with respect to the cosmic time co¬ 
ordinate t. In the homogeneous and isotropic limit (i.e. 

—>■ a^{t)5ij) we get q{t,x) —>■ —aaja} since, by defi¬ 
nition, K = Kij and TiK^ = KfKj > 0. From the 
(00) component of the Einstein equations (contracted in 
the form i?)) = 87rG[T^ — TS'l^/2]) the equation obeyed 
by the deceleration parameter is: 

q{t, x)TtK^ = SttGUp + 3p)/2 + 3(^/2 -h A(t, f)], 
where p and p are, respectively, the energy density and 


the pressure of the fluid. To lowest order in v?, A{t, x) is 
given by [{p -\- p)v? + K(^v? + 2r]u’^u^Kki\ where K^t = 
Kki — denotes the traceless part of the extrinsic 

curvature Q. Thanks to the momentum constraint Q 
the peculiar velocity is proportional to a combination of 
gradients of the extrinsic curvature and of its trace. We 
conclude that to lowest order in the spatial gradients the 
sign of the deceleration parameter is solely determined 
by the bulk viscosity while the shear viscosity only ap¬ 
pears to higher order and it affects physical scales that are 
much smaller than K~^ where is, up to an irrelevant 
numerical factor, the inhomogeneous generalization of 
the Hubble radius. The nonperturbative results confirm 
that the sign of q{t,x) is determined, to leading order, 
by the bulk viscosity coefficient Q. Over smaller scales 
the velocity field can be averaged in terms of a consistent 
coarsening that respects both the Hamiltonian and the 
momentum constraints. In summary, if the strong energy 
condition is enforced (e.g. p = 0 and p > 0) whenever 
the bulk viscosity coefficient vanishes (i.e. C = 0) the 
deceleration parameter is positive semidefinite over the 
whole Hubble patch and to leading order in the gradient 
expansion. Thus, from the range of validity of the ap¬ 
proximation scheme employed in [ij and from the fully 
nonlinear analysis based on the gradient expansion, we 
conclude that the shear viscosity can only affect typi¬ 
cal scales smaller than the Hubble radius. The author 
wishes to thank S. Floerchinger and U. Wiedemann for 
useful exchanges of ideas. 
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